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Abstract 

Managing interference in a network of macrocells underlaid with femtocells presents an important, 
yet challenging problem. A majority of spatial frequency/time reuse based approaches partition the users 
by coloring the interference graph, which is shown to be suboptimal. Some spatial time reuse based 
approaches schedule the maximal independent sets (MISs) of users in a cyclic, (weighted) round-robin 
fashion, which is inefficient for delay-sensitive applications. Our proposed policies schedule the MISs 
in a non-cyclic fashion, which aim to optimize any given network performance criterion for delay- 
sensitive applications while fulfilling minimum throughput requirements of the users. Importantly, we 
do not take the interference graph as given as in existing works; we propose an optimal construction 
of interference graph. We prove that under certain conditions, the proposed policy achieves the optimal 
network performance. For large networks, we propose a low-complexity algorithm for computing the 
proposed policy, and prove that under certain conditions, the policy has a competitive ratio (with respect 
to the optimal network performance) that is independent of the network size. The proposed policy can 
be implemented in a decentralized manner by the users. Compared to the existing policies, our proposed 
policies can achieve improvement of up to 130 % in large-scale deployments. 


I. Introduction 

As more and more devices are connecting to cellular networks, the demand for wireless 
spectrum is exploding. Dealing with this increased demand is especially difficult because most 
of the traffic comes from bandwidth-intensive and delay-sensitive applications such as multi- 
media streaming, video surveillance, video conferencing, gaming etc. These demands make it 
increasingly challenging for the cellular operators to provide sufficient quality of service (QoS). 
Dense deployment of distributed low-cost femtocells (or small cells in general, such as microcells 
and picocells) has been viewed as one of the most promising solutions for enhancing access to 
the radio spectrum fl] [2j. Femtocells are attractive because they can both extend the service 
coverage and boost the network capacity by shortening the access distance (cell splitting gain) and 
offloading traffic from the cellular network (offloading gain). However, in a closed access network 
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when only registered mobile users can connect to the femtocell base station, dense deployment of 
femtocells operating in the same frequency band leads to strong co-tier interference. In addition, 
the macrocell users usually operate in the same frequency, the problem of interference (to both 
femtocells and macrocells) is further exacerbated due to cross-tier interference across macrocells 
and femtocells. In this work, we study a closed access network. Hence, it is crucial to design 
interference management policies to deal with both co-tier and cross-tier interference. 

Interference management policies specify the transmission schedule and transmit power levels 
of femtocell user equipments FUEs (femtocell base stations FBSs) and macrocell user equipments 
MUEs (macrocell base stations MBSs) in uplink (downlink) transmissions. [[] An efficient inter¬ 
ference management policy should fulfill the following important requirements (we will discuss 
in Section [TT[ state-of-the-art policies do not fulfill one or more of the following requiremenst): 

• Interference management based on network topology. Policies must take into account that 
uplink transmissions from neighboring UEs create strong mutual interference, but must also 
recognize and take advantage of the fact that non-neighboring UEs do not. Hence, the 
network topology (i.e. locations of femtocells/macrocells) must play a crucial role. 

• Limited signaling for interference coordination: In dense, large-scale femtocell deployments, 
the UEs cannot coordinate their transmissions by sending a large amount of control signals 
across the network. Hence, effective policies should not rely on heavy signaling and/or 
message exchanges across the UEs in the network. 

• Scalability in large networks: Femtocell networks are often deployed on a large scale (e.g. 
in a city). Effective policies should scale in large networks, i.e. achieve efficient network 
performance while maintaining low computational complexity. 

• Support for delay-sensitive applications: Effective policies must support delay-sensitive ap¬ 
plications, which constitute the majority of wireless traffic. 

• Versatility in optimizing various network performance criteria: The appropriate network per¬ 
formance criterion (e.g. weighted sum throughput, max-min fairness, etc.) may be different 
for different networks. Effective policies should be able to optimize a variety of network 
performance criteria while ensuring performance guarantees for each MUE and each FUE. 

In this work, we propose a novel, systematic, and practical methodology for designing and 

'For brevity, we will focus on uplink transmissions hereafter. 
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implementing interference management policies that fulfill all of the above requirements. Our 
proposed policies aim to optimize a given network performance criterion, such as weighted 
sum throughput and max-min fairness, subject to each UE’s minimum throughput requirements. 
Our proposed policies can efficiently manage a wide range of interference. We manage strong 
interference between neighboring UEs by using time-division multiple access (TDMA) among 
them. We take advantage of weak interference between non-neighboring UEs by finding maximal 
sets of UEs that do not interfere with each other and allowing all the UEs in those sets to transmit 
at the same time. Specifically, we find the maximal independent sets (MISs) [^] of the interference 
graph, and schedule different MISs to transmit in different time slots. The scheduling of MISs 
in our proposed policy is particularly designed for delay-sensitive applications: the schedule of 
MISs across time is not cyclic (i.e. the policies do not allocate transmission times to MISs in 
a fixed (weighted) round-robin manner), but rather follows a carefully designed nonstationary 
schedule, in which the MIS to transmit is determined adaptively online. For delay-sensitive 
applications, cyclic policies are inefficient because transmission opportunities (TXOPs) earlier 
in the cycle are more valuable than TXOPs later in the cycle (earlier TXOPs enhances the 
chances of transmission before delay deadlines). The cyclic polices are unfair to UEs allocated 
to later TXOPs. 

Another distinctive feature of our work is that we do not take the interference graph as given 
as in most existing works; instead, in our work we show how to choose the interference graph 
that maximizes the network performance. Specifically, in our construction of interference graphs, 
we determine how to choose the threshold on the distance between two cells, based on which 
we determine if there is an edge between them, in order to maximize the network performance. 
Moreover, we prove that under certain conditions, the proposed policy, computed based on the 
optimal threshold, can achieve the optimal network performance (weighted sum throughput) 
within a desired small gap. Note that for large networks, in general it is computationally 
intractable to find all the MISs of the interference graph [3]. We propose efficient polynomial¬ 
time algorithms to find a subset of MISs, and prove that under certain conditions, the proposed 
policy, computed based on the constructed subset of MISs, can achieve a constant competitive 

2 A set of vertices in which no pair is connected by an edge is independent (IS) and if it is not a subset of another IS then it 
is MIS. 
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ratio (with respect to optimal weighted sum throughput) that is independent of the network size. 

Finally, we summarize the main contributions of our work: 

1. We propose interference management policies that schedule the MISs of the interference 
graph. The schedule is constructed in order to maximize the network performance criterion 
subject to minimum throughput requirements of the UEs. Moreover, the schedule adapts to the 
delay sensitivity requirements of the UEs by scheduling transmissions in a non-stationary manner. 

2. We construct the interference graph by comparing the distances between the BSs with a 
threshold (i.e. there is an edge between two cells if the distance between their BSs is smaller than 
the threshold). We develop a procedure to choose the optimal threshold such that the proposed 
scheduling of MISs leads to a high network performance. Importantly, we prove that under 
certain conditions, the proposed scheduling of MISs based on the optimal threshold achieves 
within a desired small gap of the optimal network performance (weighted sum throughput). 

3. It is intractable to find all the MISs in large networks, we propose an approximate algorithm 
that computes a subset of MISs within polynomial time. We prove that under certain conditions, 
the proposed policy based on this subset of MISs has a constant competitive ratio (with respect 
to the optimal weighted sum throughput) that is independent of the network size. 

The rest of the paper is organized as follows. In Section [TT] we discuss the related works 


and their limitations followed by the system model and problem formulation in Section III 


and [TV] respectively. The design framework and its low-complexity variant for large networks 
are discussed in Section [V] and Section |VI] respectively. In Section |VII| we use simulations to 


compare the proposed policy with state-of-the-art followed by conclusion in Section VIII 


II. Related Works 

In this section we compare our proposed policy with state-of-the-art policies. Existing in¬ 
terference management policies can be categorized in two classes: 1) policies based on power 
control, and 2) policies based on spatial time/frequency reuse. 


A. Interference Management Policies Based on Power Control 

The first and most widely-used interference management policies [|4]—[ 10J are based on power 
control. In these policies, all the UEs in the network transmit at a constant power at all time 
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(provided that the system parameters remain the same) in the entire frequency band [] When the 
cross channel gains among BSs and UEs are high, simultaneous transmissions at the same time 
and in the same frequency band will cause significant interference among cells. Such strong 
interference is common in macrocells underlaid with femtocells. For example, in G3 it is 
shown that that interference from MUEs near the FBS severely affects the uplink transmissions 
of FUEs. Also, in offices and apartments, where FBSs are installed close to each other, inter¬ 
cell interference is particularly strong. In contrast, our proposed solutions mitigate the strong 
interference by letting only a subset of non-interfering UEs to transmit at the same time. 


B. Interference Management Policies Based on Spatial Time/Frequency Reuse 

Some existing works mitigate strong interference by letting different subsets of UEs to transmit 
in different time slots (spatial time reuse) p~2|-[20] or in different frequency channels (spatial 
frequency reuse) [f2Tj|-[27]. Specifically, they partition UEs into disjoint subsets such that the 
UEs in the same subset do not interfere with each other. Given the same partition of the UEs, the 
policies based on spatial time reuse and those based on spatial frequency reuse are equivalent. 
Hence, we focus on policies based on spatial time reuse hereafter. 

Some policies based on spatial time reuse, partition the UEs based on the coloring of the 
interference graph h [12]— which is not efficient. In general, a set of UEs with the same 
color (i.e. the UEs who can transmit simultaneously) may not be maximal (See Fig. 1 a), in the 
sense that there may be UEs who do not interfere but have different colors (we will also show 


this in the motivating example in Subsection |IV-B| ). In this case, it is more efficient to also let 
those non-interfering UEs to transmit simultaneously, although they have different colors. Hence, 
the partitioning based on coloring the interference graph is not efficient, because the average 
number of active UEs (i.e. the average cardinality of the subsets of UEs with the same color) is 
low. 

Some policies based on spatial time reuse [17]-[j20]] partition the UEs based on the MISs of 
the interference graph, which is more efficient, because we cannot add any more UEs to an MIS 
without creating strong interference. However, they are still inefficient compared to our proposed 


’Although some power control policies 14] go through a transient period of adjusting the power levels before converging to 
the optimal power levels, the users maintain the constant power levels after the convergence. 
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policies for delay-sensitive applications. Specifically, they schedule different MISs in a cyclic 
and (weighted) round-robin manner, in which each UE transmits at a fixed position in each 
cycle. For delay-sensitive applications, earlier positions in the cycle are more desirable because 
they enhance the chances of transmitting prior to delay deadlines. Hence, a cyclic schedule is 
not fair to the UEs allocated to later positions. In contrast, our proposed policies schedule the 
MISs in an efficient, nonstationary manner for delay-sensitive applications. 

Another notable difference from the existing works based on spatial time/frequency reuse is 
that they usually take the interference graph as given. On the contrary, our work discusses how 
to construct the interference graph optimally such that the network performance is maximized. 


C. Other Interference Management Policies 


Besides the above two categories, there are several other related works. For instance in [281 
|, the authors propose reinforcement learning and evolutionary learning techniques for the 
femtocells to learn efficient interference management policies. In [[28]], the femtocells learn the 


fixed transmit power levels, while in [29], the femtocells leam to randomize over different 
transmit power levels. However, the interference management policies in [28] and p9| cannot 
provide minimum throughput guarantees for the UEs. In constrast, we provide rigorous minimum 


throughput guarantees for the UEs. In both [281 [29] the femtocell UEs need to limit their 
transmission powers in every time slot such that the SINR of the macrocell UE is sufficiently 
high. If there is strong interference between some femtocells and the macrocell, the femtocell 
UEs will always transmit at lower power levels, leading to a low sum throughput for them. 


Another method to mitigate interference is to deploy coordinated beam scheduling [30] [31]. 


In [30] and [31], the authors schedule a subset of beams to maximize the total reward associated 
with the scheduled subset, where the reward per beam reflects the channel quality and traffic. 
The first difference from our work is that the proposed approach schedules a fixed subset of 
beams and leaves the other UEs inactive. Hence, some UEs have no throughput, which means 
the minimum throughput as well as the delay-sensitivity of the UEs is not satisfied. Second, we 
rigorously prove that our proposed policy achieves good performance with low (polynomial-time) 
complexity , while |30] |3_1] do not. Third, the schemes in [ |30| [31] are proposed for a specific 
network performance criterion and may not be flexible enough for other network performance 


criteria (such as the sum throughput). Finally, 30] [31] do not consider delay sensitivity of the 
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Figure 1. a) Coloring based scheduling in i) schedules less than two UEs per time slot on an average, while MIS based 
scheduling in ii) is more efficient and schedules two UEs per time slot. , b) System model illustration. 


UEs. 

III. System Model 

A. Heterogeneous Network of Macrocells and Femtocells 

We consider a heterogeneous network of N femtocells (indexed by {l,2,...,iV}) and M 
macrocells (indexed by {N + 1 ,..., A r + M}) operating in the same frequency band, a common 
deployment scenario considered in practice [4j. We assume that each FBS/MBS serves only one 
FUE/MUE . Our model can be easily generalized to the setting where each BS serves multiple 
UEs, at the expense of complicated notations to denote the association among UEs and BSs. 
For notational clarity, we focus on the case where each BS serves one UE, and will demonstrate 


the applicability of our work to the setting where one BS serves multiple UEs in Section VII 


Since there is only one FUE or MUE in a femtocell or a macrocell, the index of each UE and 
that of each BS are the same as the index of the cell they belong to. We focus on the uplink 
transmissions. The proposed framework can be applied directly to the downlink scenarios in 
which each BS serves one UE at a time. See Fig. lb) for an illustration of a 3-cell network 
with N = 2 femtocells and M = 1 macrocell. Each UE i chooses its transmit power p, from a 
compact set V l C M + . We assume that 0 G Vi, Vi G (1, ...N + M}, namely a UE can choose 
not to transmit. The joint power profile of all the UEs is denoted by p = (p 3 , ....p N+M ) G V, 
where V = flf i M P> - The power profile of all the UEs other than i is denoted by p When a 
UE i chooses a transmit power pi, the signal to interference and noise ratio (SINR) experienced 
at BS i is 7*(p) = v g “p’+rr 2 ’ here !jp, is the channel gain from UE j to BS i, and erf is the 
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noise power at BS i. Since the BSs cannot cooperate to decode their messages, each BS i treats 
the interference as white noise, and gets the following throughput [4j] r i(p) = log 2 (l + 7i(p))- 

B. Interference Management Policies 

The system is time slotted at t = 0,1,2..., and the UEs are assumed to be synchronized as 


in [32]. At the beginning of time slot t, each UE i decides its transmit power p\ and obtains a 
throughput of rf p 4 ). Each UE i’s strategy, denoted by 7r* : Z + = {0,1,..} —> Vi, is a mapping 
from time t to a transmission power level />,, G V, . The interference management policy is then the 
collection of all the UEs’ strategies, denoted by n = ( 7 r l5 ..., n N+M )- Each UE is delay sensitive 


and hence discounts the future throughputs as in [331. The average discounted throughput for UE 


i is given as Rf n) = (1 — 5) ^ r i( p 4 ), where p* = (7Ti(t),..., n N+M (t)) is the power profile at 

t =0 


time t, and 5 G [0,1) is the discount factor assumed to be the same for all the UEs as in [32], 


133]. We also assume the channel gain to be fixed over the considered time horizon as in [18], 


1191, [122] [32], [33]. We will illustrate in Section VI-D and Section VII-C that the proposed 


framework can be adapted to the scenarios in which the channel conditions are time-varying. 

An interference management policy 7v const is a policy based on power control [4j], if 7 r(f) = p 
for all t. Write the joint throughput profile of all the UEs as r(p) = (r’ 1 (p), ■■■, Uv+m(p))- Then 
the set of all joint throughput profiles achievable by policies based on power control can be 
written as 7 Z const = (r(p), p G V}. As we have discussed before, our proposed policy is based 
on MISs of the interference graph. The interference graph G has M + N vertices, which are 
the M + N cells. There is an edge between two cells if their cross interference is high. We will 
describe in detail how to construct the interference graph in Section V. Given an interference 
graph, we write I G = { If ,..., lf G} } as the set of all the MISs of the interference graph. Let p 7 ? 
be a power profile in which the UEs in the MIS If transmit at their maximum power levels, 
namely p k = p'f ax = max V/. if k G If and p k = 0 otherwise. Let i p MIS i G ) = { p ; ? . .... p /; (G) } 
be the set of all such power profiles. Then n is a policy based on MIS if n(t) G V MIS ( G ) 
for all t. We denote the set of policies based on MISs by n M/s ( G ) = {77 : Z + —> V MI ‘^ G )}. 
The set of joint instantaneous throughput profiles achievable by policies based on MIS is then 
P^mis{G) _ | r (p) ; p £ w e will prove in Theorem 1 that the set of joint discounted 

throughput profiles achievable by policies based on MIS is = com{TZ MIS ^}, where 

conv{ A" } representing the convex hull of set X. 
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IV. Problem Formulation 

In this section, we formalize the interference management policy design problem, and sub¬ 
sequently give a motivating example to highlight the advantages of the proposed policy over 
existing policies in solving this problem. 

A. Policy Design Problem 

The designer of the network (e.g. the network operator) aims to design an optimal interfer¬ 
ence management policy it that fulfills each UE i’s minimum throughput requirement 
and optimizes a chosen network performance criterion W(Ri(n), ...., R n+m (it)). The network 
performance criterion W is an increasing function in each R,. For instance, W can be the 
weighted sum of all the UEs’ throughput, wf UE R^ tt) + YljLi w j IUE with 

Y^iLi wf UE + i w f WE = 1 and wf ,UE . w EUE > 0 . We emphasize that the higher-priority 
of MUEs can be reflected by setting higher weights for the MUEs (i.e. yjMUE > w EUE \/i = 
{1, ..., V}, V) = {1,..., M}), and by setting higher minimum throughput requirements for 
MUEs. The performance criterion W can also be max-min fairness (i.e., the worst UE’s through¬ 
put min,; R, t (7v)) or the proportional fairness. The policy design problem is given as follows. 

Design Problem max ff W(Ri(7t), R n+m {7t)) (1) 

s.t. Ri(TT) > R™ in , Vi e (1,..., N + M} 

The key steps and the challenges in solving the design problem are as follows: 1) How to 
determine the set of achievable throughput profiles? Note that the set depends on the discount 
factor 5. It is an open problem to determine the set of achievable throughput profiles, even for the 
special case of d = 0 (i.e. the set of throughput profiles achievable by policies based on power 
control). 2) How to construct the optimal policy that achieves the optimal target throughput 
profile? The optimal policy again depends on 5. It is much more challenging to determine the 
policy for delay-sensitive applications (i.e. <5 < 1) than for delay-insensitive applications (i.e. 
5 —» 1), because the optimal policy is not cyclic. 3) How to construct a distributed policy with 
minimum communication overhead? 

B. Motivating Example 

We consider a network of 5 femtocells. On the left plot of Fig. 1 a), we have portrayed the 
interference graph of this network. Each vertex denotes a pair of FBS and its FUE. Each edge 
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denotes strong local interference between the connected vertices (i.e. the distance between the 
FBSs is below some threshold). The interference graph is a pentagon, where each UE interferes 
only with two neighbors. We show the partitioning of the UEs by coloring the interference graph. 
There are three colors, and there is one color (i.e. black) to which only one UE belongs. On the 
right plot of Fig. 1 a), we show the 5 MIS’s, each of which consists of two UEs. Note that the 
MIS are not disjoint. For illustrative purposes, suppose that the 5 femtocells and their UEs are 
symmetric, in the sense that all the UEs have maximum transmit power of 30mW, direct channel 
gain of 1, cross channel gain of 0.25 between the neighbors, noise power at the receiver of 2mW, 
minimum throughput requirement of 1.2 bits/s/Hz, and discount factor of 0.8 representing delay 
sensitivity. For simplicity, we set the cross channel gain between non-neighbors to be 0. 

We compare our proposed policy against the following policies discussed in Section |TT[ 

. Policies based on power control [4j [34], in which each UE chooses a constant (time- 
invariant) power level all the time. 


Coloring-based TDM A policies [12] [16], in which the UEs are partitioned into mutually 
exclusive subsets by coloring the interference graph; in each time slot, all the UEs of one 
color are chosen to transmit. In this example, 3 colors are required and there will always 
exist a color to which only one UE belongs. Hence, the average number of active UEs in 
each time slot is less than 2. Note that the optimal performance of coloring based frequency 
reuse policies is the same as the optimal performance that can be attained by any coloring 
based TDMA of any arbitrary cycle length. This is due to the fact that FDM and TDM are 
equivalent provided the frequency/time can be divided arbitrarily. 


• Cyclic MIS-based TDMA policies [18] [19], in which different MISs of UEs are scheduled 
in a cyclic manner. In this example, there are 5 MISs, each of which consists of 2 UEs. 
Hence, the average number of active UEs in each time slot is 2. This is the major reason 
why MIS-based TDMA policies are more efficient than coloring-based TDMA policies. To 
completely specify the policy we must also specify a cycle length and order of transmissions; 
note that the efficiency of the policy will depend on the cycle length due to delay sensitivity. 
We illustrate the performance of the above policies vs the proposed policy in Table 1. Remarkably, 
the proposed policy is not only much more efficient than existing policies, it is much easier to 
compute. To compare with constant policies, note simply that finding the optimal constant policy 


April 28, 2015 


DRAFT 








11 


Table I 

Comparisons against spatial reuse TDMA based policies 


Policies 

Max-min throughput (bits/s/Hz) 

Performance Gain % 

Optimal constant power 

1.32 

21.2% 

Optimal Coloring TDMA (arbitrary L) 

1.33 (Upper Bound) 

20.3 % 

Optimal MIS TDMA (L=5) 

1.36 

17.6 % 

Optimal MIS TDMA (L=7) 

1.49 

7.8 % 

Optimal Proposed 

1.60 

- 


is intractable [34] in general, because the optimization problem is non-convex due to the mutual 
interference. To compare with different classes of TDMA policies, note that for (coloring-based 
and MIS-based) cyclic TDMA policies, the complexity of finding the optimal cyclic policy of 
a given length grows exponentially with the cycle length (and exponentially with the number 
of MISs when the cycle length is large enough for reasonable performance). To get a hint of 
why this is so, note that in a cyclic policy, the UE’s performance is determined not only by the 
number of TXOPs in a cycle but also by the positions of the TXOPs since UEs are discounting 
their future utilities (due to delay sensitivity). Thus, it is not only the length of the cycle that is 
important but also the ordering of transmissions within each cycle. For instance, for the 5-UE 
case above, achieving performance within 10% of the optimal nonstationary policy requires that 
the cycle length L be at least 7, and so requires searching among the thousands (16800) |^] of 
different nontrivial schedules (the schedules in which each UE transmits at least once in each 
cycle) of cycle length 7. Even this small problem is computationally intensive. For a moderate 
number of 10 femtocells, assuming a completely connected interference graph which has 10 
MISs, and a cycle length of 20, we need to search more than ten billion (i.e. 10 10 ) non-trivial 
schedules - a completely intractable problem. 


V. Design Framework 

In this section, we develop the general design framework for solving the design problem 
(1). We will provide sufficient conditions under which our proposed framework is optimal, and 

4 We compute the number of nonstrivial schedules by exhaustively searching among all the possible policies. 
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Figure 2. Steps in the design framework. 


demonstrate a wide variety of networks that fulfill the sufficient conditions. 


A. Description of the Proposed Design Framework 

The proposed methodology for solving the design problem consists of 5 steps which are 
illustrated in Fig. [2] . We describe them in detail as follows. 

1) Step 1. The Designer Gathers Network Information: The designer is informed by each 
BS i of the minimum throughput requirement Rf m of its UE, the channel gain from each UE 
j to its receiver g Jt , its UE’s maximum transmit power level pf ax , the noise power level at its 


receiver of, and its location as in [22] [231. Such information is sent to the designer via the 
backhaul link. In some circumstances, the information about the location of FBSs is available 
to the femtocell gateways [23J, who can send this information to the designer. 

2) Step 2. The Designer Constructs the Interference Graph and Computes the MISs: The 
designer constructs the interference graph using the information of cell locations obtained in Step 


1. Specifically, it uses a distance based threshold rule as in [ 131 [35] to construct the graph: 
there is an edge between two cells if the distance between BSs in these two cells is smaller than 
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a threshold D [^] . Given the threshold D, we denote the resulting graph by Go, and the set of 
its MISs by 1 Gd , which can be calculated as in [3]. We assume that the distance threshold D is 
fixed for now, and will discuss how to select the threshold in the next subsection. 

3) Step 3. The Designer Characterizes Achievable Throughput: Based on the MISs com¬ 

puted in Step 2, the designer identifies the set V mis ^ Gd \8) of throughput vectors achievable 
by MIS-based policies. Note that \J mis ( Gd ')(S) depends on the discount factor. Recall that 
j^mis(g d ) _ | r (p) : p (= pM/s(G D ) j j s set 0 f instantaneous throughput profiles achievable 
by MIS-based policies in yI mis ^ Gd \ The theorem below proves that \? mis ( g d) (<j) j s a convex 
hull of 7 ^ mi s(g d ), p e yMis(G D ) w j ien discount factor 5 > 1 — where s(Gd) is the 

number of MISs in the interference graph Go- 

Theorem 1: Given the interference graph Go, for any s(G D y the set of throughput 

profiles achieved by MIS-based policies is V M/s ^ Go ^(5) = y MIS ( G D)_ 

See the detailed proof in the Appendix at the end. 

Proof Sketch: The main step involved in proving the above is to derive the conditions on the 
discount factor such that each throughput vector in \? mis ( g d) can |-, e decomposed into a current 
throughput vector which belongs to 7 Z mis( - Gd ) and a continuation throughput which belongs to 
\>mis{g d ) ' y 0 derive the conditions, we show that for any vector in \? mis ( g d) there exists at least 
one throughput vector in 'JZ MIS( - Gn ) t 0 decompose the vector. Since the continuation throughput 
also belongs to V M/s ( Gfl ), it can be decomposed as well in a similar fashion. Hence, all the 
vectors in \? mis ( g d) are achievable. 

Theorem 1 is important, because it analytically characterizes the set of throughput profiles 
achievable by MIS-based policies, and gives us the requirements that need to be fulfilled by the 
discount factor. 

4) Step 4. The Designer Determines the Optimal Target Weights: Among all the achiev¬ 
able throughput profiles identified in Step 3, the designer selects the target throughput profile 
in order to optimize the network performance. Note that each UE Ts average throughput R, 
can be expressed as a convex combination of the instantaneous throughput vectors achieved by 
MIS-based policies (i.e. the throughput vectors in 1Z m,s{Gd) ). Thus, determining the optimal 

5 Note that the interference actually depends on the distance between a BS and a UE in another cell, instead of the distance 
between two BSs. When the distance from a BS to its UE is small, then the distance between BSs is an accurate representation 
of interference. 
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target vector and its corresponding coefficient can be formulated as the following optimization 
problem : 

max W(ui(Gd), Dn+m{Gd)) 

v,ol 

s.t. yi (G D ) > FC in , V* e {1 ,...,N + M} 

s(Gd) g 

Vi(G D ) = J2 a M P^ D ),Vi e {1, ,N + M} (2) 

3 =1 

s(G d ) 

Y aj = 1, ay > 0, Mj e {1, ..,s(G d )} 
j =i 

The above optimization problem is a convex optimization problem and is easy to solve if W 
is concave (e.g. weighted sum throughput or max-min fairness). The resulting optimal target 
vector and its corresponding coefficient is given as y *(Gd) = [?/J(Gd),....,i/^ +m (Gd)] and 

cx*(G d ) = [a^G'/j)),., G*( ( ; n) {Gn). respectively. Note that the optimal value depends on the 

interference graph Go which we assume to be fixed in this section. The optimal coefficient for 
the i th MIS if , i.e a*(Go) can be interpreted as the fraction of time for which /f 'Transmits. 

5) Step 5. Each UE Implements the Policy Distributedly to Achieve the Target: The designer 
informs each UE i of the optimal coefficients, i.e. a*(Go) and the indices of MISs that UE i 
belongs to. The designer can send the above information to each BS i, who will forward the 
information to its UE. Each UE i executes the policy in Table 2. The policy in Table 2 leads to 
a non-stationary scheduling of the MISs. Note that each UE i computes its own policy online 
without information exchange. Hence, the computed policy is implemented in a decentralized 
manner by the UEs. Next we state the condition under which the policy indeed converges to the 
target vector y *(Gd)- 

Theorem 2: For any the policy computed in Table 2 achieves the target 

throughput profile y *(G D ). 

See the detailed proof in the Appendix at the end. 

Proof Sketch: We show that when s(G D )’ the policy developed in Table 2 ensures 

that the decomposition property given in Proof Sketch of Theorem 1 is satisfied in each time 
slot. This is used to show that the distance from the target, y *{Gd) strictly decreases in each 
time slot. 
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We briefly discuss the intuition behind our proposed policy. We determine which MIS to 
transmit based on a metric that can be interpreted as the “fraction of time slots allocated to an 
MIS in the future”: the MIS that has the maximum fraction of time slots in the future, i.e. the 
highest metric, will transmit at the current time slot. The metric is updated in each time slot as 
follows: the fraction of time slots for the MIS who has just transmitted will decrease, and those 
of the other MISs will increase. 

B. Constructing Optimal Interference Graphs 

In Step 2 of the design framework, we construct the interference graph by comparing the 
distances between two BSs with a threshold D. Here we show how to choose the optimal 
threshold D* and hence the optimal interference graph Go*, based on which the proposed policy 
achieves the highest network performance achievable by any MIS based policy in yI mis ( g d)_ 
Formally, the designer chooses the optimal threshold D* that results in the optimal interference 
graph G D * = argmax Gj3e g; W(y*(G D )), where Q is the set of all possible interference graphs 
constructed based on the distance rule. The designer solves the above optimization problem by 
performing Steps 2-4 for each of the \Q\ — J interference graphs as shown in Fig. 2. and chooses 
the optimal one. Note that the number \Q\ of all such interference graphs is finite and upper 
bounded by ( M + N )-G' I+N - 1 ) _|_ ^ because the number of different distances between BSs is finite 
and upper bounded by - ai +n)-(m+n-i) _|_ \ Note that the Steps 3-5 of our design framework 
can be used for any given interference graph, which is not necessarily constructed based on the 
distance based threshold rule. We assume a distance based threshold rule as a concrete example, 
in order to describe how to choose the optimal interference graph. 

C. Optimality of the Proposed Design Framework 

Our proposed design framework first constructs the interference graph based on the distances 
between BSs, and then schedules the MISs of the constructed interference graph. Then our 
proposed policy let the UEs in the scheduled MIS to transmit at their maximum power levels. 
To some extent, the interference graph is a binary quantization of the actual interference (i.e. 
“no interference” among non-neighbors and “strong interference” among neighbors). Hence, the 
performance of the proposed policy depends crucially on how close the interference graph is to 
the actual interference pattern. If we choose a smaller threshold D, the interference graph will 
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have fewer edges, the non-neighboring UEs will have higher cross channel gains. Hence, the UEs 
in a MIS may experience high accumulative interference from the non-neighbors. If we choose 
a higher threshold D, the interference graph is more conservative and will have more edges. 
Hence, some UEs outside a MIS may cause low interference and should be scheduled together 
with the UEs in the MIS. Our proposed policy will achieve performance close to optimality, if 
the interference graph is well constructed such that: 1 ) neighbors have strong interference, and 
2) non-neighbors have weak interference. Next, we analytically quantify the above intuition and 
provide rigorous conditions for the optimality of the proposed design framework. 

Let W* denote the optimal network performance, namely the optimal value of the design 
problem (1) with the performance criterion being the weighted sum throughput. We give con¬ 
ditions under which the proposed policy can achieve within e of the optimal performance 
W*. We first quantify strong local interference among neighbors as follows. Define r'(p) = 


log 2 (l + 


9iiPi 


), where MfGjf) is the set of neighbors of i in G D and let 7Z 


const 


{u(p), p € V} , 7 ^mis(g d ) _ p e -pMis(G D )j an( j yMis(G D ) _ conv | 7 Za IIS( ' GD ' > }. Note 

that r'(p) is not the actual throughput r t (p), because we do not count the interference from 
non-neighbors in r'(p). 

Definition 1 (Strong Local Interference): The interference graph Gn exhibits Strong Local 
Interference (SLI) if dominates TZf nst , in the sense that every throughput profile in 

7 Zfnst is weakly Pareto dominated by a throughput profile in ,s(Gn \ 

Definition 1 states that for an interference graph with SLI, it is more efficient to use MIS-based 
policies than constant power control policies. Next, we quantify the weak interference among 
non-neighbors. 

Definition 2 (Weak Non-neighboring Interference): The interference graph Gd has e — 

Weak Non-neighboring Interference (e-WNI) if each UE i’s maximum interference from its 
non-neighbors is below some threshold, namely Intf ax (G o) = Ylj^j\ti(G D ) HGGf QX < (2 e — 

1K 2 , Vie TV+ M}. 


Now we state Theorem 3 which uses the above two definitions to ensure optimality. 
Theorem 3: If the constructed interference graph Go* exhibits SLI and e-WNI, then the 
proposed policy computed through Steps 1-5 of Subsection V-A achieves within e of the optimal 
network performance W*. 

See the Appendix at the end for the detailed proof. 
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Proof Sketch: The set of throughput vector achievable by any policy is conv{'R const }. Denote 
the optimal throughput vector by v* G conv{lZ const }, namely W(v*) = W*. There must exist 
a vector v G conv{TZf mst } such that v > v*, because we do not count the interference from 
non-neighbors when we calculate r'(p) G r Rff nst . SLI indicates that there exists a vector V G 
yMis(c D *) suc ^ that V > v > v*. This condition implies that if hypothetically there was zero 
interference from non-neighbors, then MIS based policies will achieve the optimal throughput 
vector. However, since there is interference from non-neighbors, we use e-WNI to bound the 
loss in throughput caused by the interference from non-neighbors. Using e-WNI we can find 


a throughput profile v G \? mis ( g d*) which is within e from v' G V, 


MIS(G d *) 


. Hence, we have 


v > v > v — e and v, > R™ m — e . This shows that we can achieve a throughput vector that is 
e close to the optimal one, i.e. v > v* — e . 

Example: Consider 3 UEs and their corresponding FBS located on 3 different floors as shown 
in Fig. 3. Each UE can transmit at a maximum power of 100 mW. The channel model for 
determining the gain from a UE i to BS j, which includes the attenuation from the floor, is 


set based on [36]. Specifically, we have G u = 0.5, G 3i = 0.25 for | j — i\ — 1, Gj t = 0.0032 
for \j — i\ =2, and the noise power of 2 mW. We aim to maximize the average throughput 
while fulfilling a minimum throughput requirement of 1.2 bits/s/Hz for each FUE. Under three 
different thresholds D, we have the following three interference graphs (there are only three 
interference graphs because there are only three different values of distance between the BSs): 
1) the triangle graph {D > 4m}, 2) the chain graph (2m < I) < 4m} and 3) the edge-free 
graph {0m < D < 2m}. For each of these graphs, we apply the design framework described in 
Subsection V-A to obtain the corresponding policy, and achieve the following average throughput: 
1) 1.56 bits/s/Hz 2) 2.7 bits/s/Hz and 3) 1.5 bits/s/Hz. Hence, the chain graph is the optimal 
choice among the three graphs. Also the chain graph exhibits SLI as illustrated in Fig. 3. and 
also exhibits e— WNI for e = 0.2. Hence, the proposed policy calculated based on the chain 
graph yields an average throughput within e = 0.2 of the optimal solution W* to the design 
problem in (1) (i.e. W* < 2.9 bits/s/Hz). 


D. Complexity for computing the policy 

We only compare the computational complexity of the proposed policies against cyclic MIS- 
based TDMA policies, since determining the optimal constant power based policy is a non-convex 
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FUEL'S throughput 



Figure 3. An example to illustrate the optimality of proposed framework. 


problem and has been shown to be NP-hard [34]. We compare the two for a given interference 
graph Gd■ Both the optimal cyclic MIS TDMA and the proposed policy need to compute 
the set of MISs. Determining all the MISs in general requires large computational cost [3|. 
However, the computational complexity is acceptable if the network is small, or if the number 
of MIS, s(Gd) = 0((N + M) c ), c > 1 is bounded by a polynomial function in the number of 
vertices in G D . We will develop an approximate algorithm to compute only a subset of MISs 


within polynomial time and with performance guarantees in Section VI In our framework, the 
remaining amount of computation (other than computing MISs) is dominated by the amount of 
computation performed in Step 4, because in Step 5, the policy is computed online with a small 
amount 0{s{Gd)) of computations per time slot. In Step 4, we solve the optimization problem 
in (5) with the objective function W and linear constraints. When W is linear (e.g. weighted 
sum throughput) or is the minimum throughput of any UE (in which case the problem can be 
transformed into a linear programming), the worst-case computational complexity for solving 


(5) is 0((s(Gd) + N + M) 35 B 2 ) [37] where B is the number of bits to encode a variable. 
In contrast, the complexity of computing the optimal cyclic MIS-based TDMA policy of cycle 
length L scales by [s{Gf))] L . The complexity quickly becomes intractable when cycle lengths 
are moderately higher than N + M, which is usually needed for acceptable performance. In 
summary, the complexity of computing our policies is much lower than that of computing cyclic 
MIS-based TDMA policies. 
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E. Impact of the density of femtocells and macrocells 

The density of the network is defined as the average number of neighbors of a UE in the 
interference graph. To obtain sharp analytical results, we restrict our attention to a class of 
interference graphs with N + M vertices and H cliques of the same size. Note that a clique is 
a subset of vertices, where any two vertices are connected. Assuming that no two cliques are 
connected, we can compute the density as d = N ~^ M — 1. When the total number N + M of UEs 
remains the same and the density d increases, the number H of cliques will decrease. Since the 
vertices in a MIS can only come from different cliques, the number of MISs decreases as H 
decreases. As a result, the complexity of the policy will decrease. When the density increases, 
the multi-user interference increases, leading to a decrease in the throughput and in the network 
performance. 


VI. Efficient Interference Management With Provable Performance 
Guarantees for Large-Scale Networks 

A. Efficient Computation of A Subset of MISs 

In our design framework proposed in Section [Vj we require the designer to compute all the 
MISs in Step 2. However, computing all the MISs is in general computationally prohibitive 
for large networks. We propose an approximate algorithm to compute a subset of MISs for a 
given interference graph Go in polynomial time and provide performance guarantees for our 
algorithm. Note that the graph Go belongs to the class of unit-disk graphs [38]. 

The subset of MISs are computed as follows. 

i). Approximate Vertex Coloring: The designer first colors the vertices |^] of interference graph 


Go using the approximate minimum vertex coloring scheme in [38]. Let C\ = {1,.... C(G r>)} 
be the indices of the colors. It is proven in 


that the number of colors used is bounded by 
C*(Gd) < C{Gd) < 3 C*{Gd) where C*(Gd) is the minimum number of colors that can be 
used to color the vertices of Go- 

ii). Generating MISs in a Greedy Manner: Each color i corresponds to an independent set 
/•. For each independent set /■, the designer adds vertices in a greedy fashion until the set is 


6 


In minimum vertex coloring the objective is to use minimum number of colors and each vertex has to be assigned atleast 


one color and no two neighbors are assigned the same color. 
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maximally independent. The procedure is described in Table 3. Let the output MIS obtained 
from Table 3 be I'ffy where k(i) is the index of the MIS in the original set of MISs I 6, \ We 
write all the MISs augmented from l[, i — 1,..., C(Gd) as ■■■Ik{C(G D ))}- 

iii). Generating the Approximate Maximum Weighted MIS : Define a weight corresponding to 
each UE/vertex i as w t = r"' ox , where r r [ MX is the maximum throughput achievable by UE i 
when all the other UEs do not transmit. Given these weights, the designer ideally will like to find 
the maximum weighted MIS, namely the MIS with the maximum sum weight of its vertices. 


However, finding the maximum weighted MIS is NP-hard [39]. Hence, the designer will find 
the /^-approximate maximum weighted MIS, denoted Ik(C(G D )+iy us ’ n 8 the algorithm in [40]. 

The set of MISs computed from the above steps is then lap prox = {^kti)^ ■■■^k(C(G D )+i)}- 
Note that ■■■H(C(G D ))} ensure that all the UEs are included in the scheduled MISs, and 

^k(C(Gn)+\) ’ s included for performance improvement. Given this subset of MISs, we can define 

jMIS{G d ) 


j G D 


t G D 


Vapprox D> = [p Jfe W, ...,p Jfe(C(G o )+1) ], 7 lapprox D) = (r(p), P G Vapprox D> } and Vapprox 

com{JZ^,rox D) }- Let U approx (G D ) = {n : Z + ->■ Vapprox d) } be the set of policies in which 
only the subset of MISs are scheduled. Steps 3,4 and 5 of the design framework in Section 
|V| are performed given this subset (See Fig. 2). The results of Theorem 1 and 2 still apply 
to the policies in Tl approx (G D ) and the set of achievable throughput profiles is Vapprox ^ given 
the S > 1 — C (q d } +1 (See Corollary 1 and 2 in Appendix at the end). The target vector in 
Vapprox ^ and the corresponding coefficient is computed as in Step 4 of Section [v] and is 
denoted as y * approx (G D ), a* approx {G D ) respectively. The coefficient vector ct* approx (G D ) along 


approx \ 

with the indices of the MISs that UE i belongs to is transmitted to the BS i as in the Step 5 of 
Section [V] 

The main intuition for the procedure developed above is as follows. Steps i) and ii) find 
MISs that contain all the UEs, and hence ensure that the minimum throughput requirements are 
satisfied. Step iii) finds the MIS that contains UEs with higher weights to optimize performance. 
Given the MISs obtained in Steps i)-iii) the Steps 3-5 of the design framework are performed. 

B. Performance Guarantees for Large Networks 

In this subsection, we consider the network performance criterion as the weighted sum through¬ 
put, and give performance guarantees for the policy when we compute the subset of MISs by 


Steps i)-iii) in the Subsection VI-A Note that as we will show in the Section VIE the subset of 
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Require: Target weights a*(G) = [oJ(Gd), a *^Q D - ) (Gd)] 


Require: V — {1,.., N + M} set of vertices, w 

Initialization: Sets t = 0, cr, = a* for all j g (1,s(Gc)}. 


vector of weights of vertices, Independent set I i and 

repeat 


Adj(7 i ) where Adj(X) is the set of neighbors of X 

Finds the MIS with the maximum weight 


Initialization: = l' i , = V n (I- U Adj(/))) c , 

r* = argmax je{li .. j3(GD)} atj 


here (X) c is the complement of X 

if i g then 


While( A f i ^ 0) 

Transmits at power level p\ = pV 1 ^ 


Jf i = sort (AT ) , sort the vertices 

end if 


in Jf i in the decreasing order of the weights Wj 

Updates aj for all j G {1, s(Gd)} as follows 


v = M i ^ , here Jf i is the first vertex in 

a,. = 


I g d _ jG d u ' 

1 k(i) _ 1 k(i) u v 

^ Vj + r * 


v( =N- n ({>/} U Adj({n'}) c 

t 4 — t 1 


end 

until 0 




Table II Table III 

Algorithm run by each UE Algorithm run by the designer 


MISs perform well in large networks for other network performance metrics as well. In particular 
the performance guarantee implies that the performance scales with the optimum W(j*{Gd)) 
as the network size N + M increases. Define D EE as the distance from LIE-/' to BS-j. We make 
the following homogeneity assumption, p r . nax = p max ^ a 2 — 0 -_ R r / lvn = R min , max, D EE < A 
and Wi = v 1 u Q Here A is fixed and does not depend on the size of the network. We fix these 
parameters in order to understand the performance guarantee as a function of the network size. 
Let the channel gain g t j = rjjrrpip, where rip is the path loss coefficient. 

pinax 

We choose the trade-off variables p, (, k that satisfy p+1 < min{— 2 3i ffn 2Ctr2 , pf+R, l°g 2 (l + 
&np* 2 )} and 0 < k < 1. Any eligible triplet p, £, k will define a class of interference graphs 
that exhibit £-WNI and have maximum degrees upper bounded by p. Note that such interference 
graphs can have arbitrarily large sizes (see the example at the end of this subsection). Then the 
following theorem provides performance guarantees for the policy described in Subsection VI-A 
for this class of interference graphs. 

Theorem 4: For any interference graph that has a maximum degree no larger than p and 

7 We can extend our result to a heterogeneous network with p > p max , a 2 < a 2 , R™ ln < maxi < A and 

Wi 'T N 1 M with c as a constant. But we do not show this general result to avoid overly complicated notations. 
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pff ICIX 

exhibits C-WNI with p + 1 < min{ } , 77777 ) log 2 (l + £^ 2 )} the policy in Sub¬ 
achieves a performance W(y* (G D )) with a guarantee that W(y* (G fl )) > 


section 


VI-A 


(1 ( l+i) K) • W(y*(G D )), where 7 = ( 3 (p + 1))^^™^. 


p" 

A ri P2C,r 2 y 


See the Appendix at the end for the detailed proof. 

pinax 

Proof Sketch: The condition that the graph does not have a degree more than p < — 3 3 i ffn 2Ca2 - 


land the £-WNI condition ensure that the algorithm proposed in Subsection VI-A yields a 
feasible solution satisfying each UE’s minimum throughput constraint. Also, it is shown that 
the minimum coefficient/fraction of time allocated to ^(C{g d )+i)^ s a a PP rox,k(C(G D )+ 1 ) > (1 — 7). 
Then it is shown that if UEs in I^c(g d )+ 1 ) were t0 transmit all the time then the competitive 
ratio achieved is no smaller than 777 ) • This combined with minimum coefficient of I^c(G D )+i) 
leads to the competitive ratio guarantee of no less than ^ 77777 ^ • 

The trade-off variables p, (, k as their name suggests provide trade-offs between how large 
is the class of interference graphs for which we can provide performance guarantees, and how 
good are the competitive ratio guarantees. On one hand, a higher n allows higher p, and higher 
p and ( allow a larger class of graphs. On the other hand, as we can see from Theorem 4, higher 
p and provided that they are eligible (higher (' decrease the maximum eligible p), result in 
higher 7, and higher 7 and k give lower competitive ratio guarantees. Hence, we can tune the 
design parameters to provide different levels of competitive ratio guarantees for different classes 
of interference graphs. Next, we give an example to illustrate Theorem 4. 

Example: Consider a layout of FBSs in a K x K square grid, i.e. K 2 FBSs with a distance of 
5m between the nearest FBSs, and assume that each FUE is located vertically below its FBS at a 
distance of 1 m. Fix the parameters p max = 100 mW, a 2 = 3 mW, R mm = O.lbits/s/Hz, 77 = 0.1, 
np = 4 and the threshold D — 7 m, which gives us the upper bound p = 4 on the maximum 
degrees. We can also verify that the interference graphs under any number K 2 of FBSs exhibit 
C-WNI with ( = 0.15. Given p = 4 and ( = 0.15, we choose the minimum k = 0.17, which 
provides the highest competitive ratio guarantee of 0.53. This performance guarantee holds for 
any interference graph of any size K. 


To construct an efficient interference graph similar to Subsection V-B when the number of 


MISs are large, the designer must use the low complexity framework in Subsection VI-A In this 


case the designer computes the subset of MISs as described in Subsection VI-A and compares the 
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optimal solution obtained to decide the best threshold . Formally stated, the designer computes 


G 


= arg rnaxceg W{y 


approx 


(G)). See Fig. 2. for a comparison of the design framework 


in Subsection VI-A for large networks with that in Subsection V-A for small networks. 


C. Complexity for computing the Subset of MISs 

We show that the proposed approximation method for computing the subset of MISs described 
in Subsection |VI-A has a complexity bounded by a polynomial in the number of vertices, i.e., 


0((N + M) C ), c > 1. This is because Steps i). and iii) use the algorithms developed in [38] and 


for which the complexity has been proven to be polynomial and ii). uses a greedy strategy 
in which there can be a maximum of N + M iterations since atleast one vertex is always removed 
from J\f- in each iteration. The worst possible number of computations in an iteration is bounded 
by (N + M) 2 . Hence the upper bound of the complexity of Step ii). 0((N + M) 3 ). Hence, the 
subsets of the MISs can be computed within polynomial time, and the policy computed using 


this subset can guarantee a constant competitive ratio as shown in Subsection VI-B 


D. Extensions 

1) Construction of interference graphs based on other rules: Our design frameworks in 
Section V-A and Section VI-A do not rely on a specific method for constructing the interference 
graph. In Step 2 of the design frameworks (i.e., the step in which the interference graph is 
constructed), we can replace our distance-based construction of the interference graph with 


construction based on other criteria, such as SINR, interference levels [22], etc. Then we can 
use the resulting interference graph as the input to Step 3. For construction rules based on other 
criteria, we can also use the procedure described in Section V-B to optimize the construction 
rule (e.g., to choose the optimal threshold of SINR or the interference level, above which an 
edge is drawn between two nodes). 

Note that in the design framework in Section VI, we find a subset of MISs, instead of all MISs, 


because the network is large. To find this subset, we use the coloring algorithm in [40], which 
is known to have polynomial-time complexity for unit-disk graphs. This is where we used the 
fact that the interference graph is constructed based on distances (such that the resulting graph 
is a unit-disk graph). However, we can use other polynomial-time coloring algorithms if the 
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HI 


1.25 m 


Length = 5 m 


Fig. 4 a) 


Fig. 4 b) 


Figure 4. a) Different interference graphs for the 3 x 3 BS grid, b) Illustration of setup with 3 rooms. 


interference graph is generated based on other criteria. We can use a standard greedy coloring 
algorithm as in gg. In the next step we extend the ISs obtained by coloring to MISs. We can 
do this based on Step ii) in Section VI-A. The target weights and the corresponding schedule for 
these MISs can be generated based on Section VI-A. Results about the performance guarantees 
in terms of competitive ratio (See Theorem 4) can also be extended to this case. 

2) Incorporating uncertainty in channel gains: Our design frameworks in Section V and 
Section VI can be extended to the deployment scenarios in which the channel gains are not static. 
For fast fading, we can replace the instantaneous throughput with the expected instantaneous 
throughput in our design frameworks. For slow fading, we can track the fading by regularly 
re-computing the policy. Re-computing the entire policy every time may be costly. In Section 
VII-C we show that the designer does not need to re-compute the entire policy to get considerable 
gains compared to the state-of-the-art. Specifically, the designer fixes the interference graph that 
is selected in the beginning, and only re-computes the target weights rather than re-compute 
the optimal interference graph and the corresponding target weights. We also show that the 
performance loss incurred with respect to the latter approach, which is based on an entire re¬ 
computation is limited (8%). 

3) Incorporating beamforming: We focus on the case where each UE has one antenna. 
When UEs have multiple antennas, we can easily incorporate beamforming in our framework. 
Beamforming mitigates the interference among the UEs served by the same BS. Hence, we can 
remove the edges between UEs in the same cell from the interference graph. Then we can use 
the new interference graph as the input to Step 3 of our design framework. 

VII. Illustrative Results 

In this section, we show via simulations that our proposed policy significantly outperforms ex¬ 
isting interference management policies under different performance criteria. These performance 
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Index of the graph 

a). Optimal interference graphs for different levels of cross 
interference. 



Distance between the neighboring base stations on the grid (d meters) 
b) Performance gain of the proposed policy over other policies 


Figure 5. a) Optimal interference graph selection for 3 x 3 grid, b) Performance comparison of the proposed policy for different 
size of the grid. 


gains are obtained under varying interference levels for both small and large networks. We also 
evaluate the proposed policy when the channel conditions are time-varying due to fading. In 
this case, the designer ideally needs to re-compute the optimal interference graph each time the 
channels change at the cost of a higher complexity. We show the robustness of the proposed 
policy when we choose a fixed interference graph regardless of the time-varying fading. 

In each setting, we compare with the state-of-the-art policies described in Section |D] namely 
the constant power control based policies and the cyclic MIS TDM A based policies. We do 
not compare with coloring based TDMA/Frequency reuse policies as it was already shown 


in Subsection IV-B that the MIS based TDMA policies will always lead to better network 
performance. Throughout this section, we will set the discount factor as the minimum one 
required when we use our original design framework in Section V (namely 5 = 1 — — 
according to Theorems 1-2), and the minimum one required when we use the approximate 
design framework for large networks in Section VI (namely 5 = 1 — , C / G ^ +1 ). In this way, we 
evaluate the performance of our proposed policies under the most delay-sensitive applications. 


A. Performance Gains Under Varying Interference Levels 

Consider a 3x3 square grid of 9 BSs and corresponding UEs with the minimum distance 
between any two BSs given as d (see Fig. 4 a.). Each UE i has 5 = 0.89 and a maximum 
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power of 200 mW. Assume that the UEs an the BSs are in two horizontal hyperplanes, and 
each BS is vertically above its UE with a distance of 3.16m. Then the distance from UE i to 


another BS j is DV. E = ^3.16 2 + ( Df, p) 2 , where Dfj S is the distance between BSs i and j. 
The channel gain from UE i to BS j is g i3 = . The performance criterion is the max-min 

fairness. Under different thresholds D chosen by the designer, there are 5 possible topologies of 
the interference graph, as shown in the Fig. 4 a. For each grid size d, the optimal solution to (2) 


is computed for each interference graph as described in Subsection V-B Fig. 5 a), shows that 


under different grid sizes (i.e. different interference levels), the optimal interference graph (i.e. 
the optimal threshold) changes. As the interference level increases, the corresponding optimal 
interference graph has more edges. Fig. 5 b). compare the performances of different policies 
under different grid sizes d (i.e. different interference levels), where the interference graph is the 
optimal one. We can see that the proposed policy achieves up to 67% performance gain over 
the second best policy. Through the above results, we see that 1) it is important to construct 
different interference graphs based on the interference level, and 2) the proposed non-stationary 
schedule of MISs outperform the cyclic schedules. 


B. Performance Scaling in Large Networks 

We study a dense deployment scenario to evaluate the performance gain of our proposed 
scheme over the state-of-the-art. We allow more than one UE to transmit to a single BS, and 
will increase the number of UEs associated with a BS. Consider the uplink of a femtocell network 
in a building with 12 rooms adjacent to each other. Fig. 4 b). illustrates 3 of the 12 rooms with 
3 UEs in each room. For simplicity, we consider a 2-dimensional geometry, in which the rooms 
and the FUEs are located on a line. Each room has a length of 6 meters. In each room, there 
are P uniformly spaced FUEs, and one FBS installed on the left wall of the room at a height 
of 2m. The distance from the left wall to the first FUE, as well as the distance between two 
adjacent FUEs in a room, is (| meters. Based on the path loss model in [36], the channel 
gain from each FBS i to a FUE j is 


(DVf) 2 A n H 


, where A = 10°' 25 is the coefficient representing 


the loss from the wall, and n l3 is the number of walls between FUE i and FBS j. Each UE 
has a maximum transmit power level of 1000 mW and a minimum throughput requirement of 
Pf" n = 0.05 bits/s/Hz. For each P, the designer chooses the optimal threshold to construct the 
optimal interference graph. Note that the UEs in the same room accessing the same BSs are all 
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a) Comparing the proposed policy using sum throughput as the 
metric. 



b) Comparing the proposed policy using minimum throughput 
as the metric 


Figure 6. Comparing the proposed policy against others for a) sum throughput as the metric, b) minimum throughput as the 
metric. 


connected to each other in the interference graph, since the distances between their receiving 
BSs is 0. 

We vary the number P of FUEs in each room from 5 to 15. We fix the 5 = 0.97 , i.e. the 
least value it can take based on the largest number of UEs per room, i.e. P = 15. For each 


P, the designer constructs the optimal graph G as described in Subsection VI-A using the low 
complexity method as the number of UEs is large. Under all considered values of P, the optimal 
interference graph connects all the UEs in adjacent rooms with edges and does not connect the 
UEs in non-adjacent rooms. We use the same optimal graph to compute the optimal cyclic MIS 
TDMA of cycle length L. The cycle length is varied from 12 to 58 depending upon the number 
of UEs (we try to choose as large cycle lengths as possible to maximize performance within a 
feasible computational complexity). The number of non-trivial cyclic policies under different P 
may vary from 10 8 to even more than 10 50 which renders exhaustive search to be intractable. 
Hence, for each P we do a randomized search in 4 million policies to search for the optimal one. 
Fig. 6. compares the performance of different policies in terms of both the max-min fairness and 
the sum throughput. The constant power policy cannot satisfy the feasibility conditions for any 
number of UEs in each room. The performance gain over cyclic MIS TDMA policies increases 
as the network becomes larger. When there are 15 UEs in each room, we can improve the worst 
UE’s throughput by 131% compared to cyclic MIS TDMA policies. 
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Figure 7. Illustrating the robustness of interference graph selection based on path loss. 


C. Performance under Dynamic Channel Conditions 
We consider a 9-cell network with a grid size d = 4.74 m, where each BS is vertically above 


its UE at a distance of 3.16 m as in Subsection VII-A Each UE has a maximum power level 
of 1000 mW, noise power of 1 mW and 5 = 0.89. The channel gain is the product of path loss 
as in Subsection VII-A and fading component, ~ Rayleigh(/3) . Here, we assume that the 
fading component changes every 50 time slots independently and the new channel conditions 
are reported to the designer by each FBS as in Step-1 in Subsection |V-A| The designer has the 
choice of re-computing the optimal interference graph and thereby the optimal target every 50 
time slots at the cost of a higher complexity, or choosing a fixed optimal interference graph 
based on the channel gains computed from the path loss model (which will be graph 3 in Fig. 4 
a) and selecting the optimal target every 50 time slots based on it. In Fig. 7 we compare the loss 
due to choosing a fixed interference graph with choosing the optimal interference graph every 
50 time slots. We average the performance for a duration over a total of 10000 time slots for 
a fixed f3. In Fig. 7, we see that for a low f3, i.e. /3 — 0.1 which implies a lower variance in 
fading the loss is only 1% and even when /3 is large, i.e. /3 = 1 then as well the loss is 8%. We 
also compare with Cyclic MIS TDMA , cycle length L = 9 and optimal constant power policy, 
the performance gain with the proposed policy using a fixed interference graph is consistently 
10 % for varying fading conditions /3, while choosing the optimal interference graph can lead 
to a maximum gain of 20%. 
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VIII. Conclusion 

In this work, we proposed a novel and systematic method for designing efficient interfer¬ 
ence management policies in a network of macrocell underlaid with femtocells. The proposed 
framework relies on constructing optimal interference graphs and optimally scheduling the MIS 
of the constructed graph to maximize the network performance given the minimum throughput 
requirements. Importantly, the proposed policy is non-stationary and can address the require¬ 
ments of delay sensitive users. We prove the optimality of the proposed policies under various 
deployment scenarios. The proposed policy can be implemented in a decentralized manner with 
low overhead of information exchange between BSs and UEs. For large networks, we develop 
a low-complexity design framework that is provably efficient. Our proposed policies achieve 
significant (up to 130 %) performance improvement over existing policies, especially for dense 
and large-scale deployments of femtocells. 


IX. Appendix 

We would begin by defining a self-generating set similar to what is defined in repeated 


game theory [42]. However, our definition is less restrictive since it does not involve incentive 
compatibility as required in strategic user setting in repeated games. 

Definition 1. Decomposability : A throughput vector v E M ;¥+m is decomposable on a 
set W with respect to a discount factor S, if there exists a power profile p E p MIS ( G w an d a 
mapping 7 : r p MIS ( G o) yy suc h that Vi E {1, ...N + M} 


Vi = (1 - 5)a(p) + <57*(P) 


(3) 


Define X>(W, 5) = (v : v is pure action decomposable on W}. 

Definition 2. Self-Generation : W is self-generating with respect to discount factor S, if each 
throughput vector v G W is decomposable on W, thus W C V(W, 5) 

Theorem 1: Given the interference graph Gd, for any the set of throughput 

profiles achieved by MIS-based policies is \> mis ^ Gd \8) = y A//s, ( G o)_ 

Proof: v mis ( g d) — com{lZ MIS ( GD ' ) }, where 7 ^mis{g d ) _ |r(p^' Vj E {1,..., s(Gd)}} 
and convjX} is the convex hull of the set X. We will show that is self-generating for 

5 > 8. If we can show that the set is self-generating then we can contruct an MIS-based policy 
which achieves any given target vector in that set. This is explained as follows. Let us assume 
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that V MS ( Gfl ) is self-generating with respect to certain 5, then given a vector v £ \>mis(g d ) 

tG" /7 jG j~y jG j~y 

it can be decomposed as, v = (1 — 5) r(p j ) + ^7(p J ). The vector 7 (pv ) obtained on 
decomposition is treated as the target vector for the transmissions starting the next period. We 
know that 'y(p Ij °) £ \? mis ( g d) an d yMis(G D ) j s assumec i t Q be self-generating, which means 
7 (pv ) can also be decomposed in the same manner by a certain MIS based power profile 
p 1 " £ 7 ^mis(g d )' -pbis slC p can be recursively followed for all the future periods to generate 
a policy and hence, the target v is achieved. Next, we show the conditions under which the set 
1 jmis{g d ) j s j nc i ee( j self-generating. 

Let v £ V MIS ^ and we can express v = «fc r (p /fe °)i ot k > 0, V/c £ {1, ...^s{G D )} and 

Y^k=\^ a k = 1- If V hfIS!Cn ' 1 is self-generating with respect to certain discount factor 5 then 
3 p I i ° and -y(p Ij ) £ \> mis (°d) a bout which v can be decomposed as follows 

v = (1 - <5)r(p 7 j D ) + Sj(p I J °) (4) 

Next, we come up with sufficient conditions for 7 (p 7 > " j £ \? mis ( g d) 

7(plfD) = v-a-|Mp 

E Ol k , t g d. Qtj — (1 — (5) . t g d . 

yr(p4 ) + ^—^- Mp 1 ,' ) 

k=l ,k^j 

Observe that J2l == 1. If 0 < ^ < 1, V/e £ {1,..., s(G D )}, k ^ j and 
0 < Qj ^ ^ < 1 then 7 (p 7 J °) £ \? mis ( g d) - This can be combined into one condition on <5 
given as 

5 > {max{a fc }, 1 — ay } 

t G D 

Note that decomposition condition requires for the existence of atleast one profile p j . This 
means we can choose the least possible bound on 5, which is sufficient to ensure that there will 

t G D 

exist at least one profile p j for decomposition. 

8 > . min {max{a fe }, 1 - ay}, 

je{l,..s(Go)} k^j 

= 1 - aq s (G 0 )], 

where q W g d )] = max je{ i r .. iS(GD)} {ai,...,a s ( Gfl )}. Also, Y£!i )(X k = 1> °k > 0, VA; £ 

(1,..., s(G*d)} yields that aq s (G D )] > s (c i:) y Hence, the condition <5 > 1 — is sufficient 
to ensure decomposition of every vector in the convex hull. Thus, y MIS ( G o) i s self-generating 
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for all the discount factors 5 > <5 = 1-Hence, we have been able to show that for 

a discount factor 5 > 1 — each vector in \? mis ( g d) can achieved by an MIS based 
policy, i.e. \> mis ( g d) c V mis ^ Gd ^(S). The throughput achieved by any MIS based policy is 
R( 7r) = (1 — 5) J2^L 0 8 t r(p t ), here r (p t ) G 'JZ mis ( g d)_ Since the coefficients (1 — 5)5* > 0 and 
the sum of the coefficients of the throughput vector sum to 1 , i.e. — = 1 > this implies 

R(n) G yMis(G D ) Hence, V MIS ^ G °\8) C y MIS (G D )_ Therefore, \j> mis ( g d) (j) — yMis(G D ) f or 

5 - 1 _ rtko)’ Next, we give a corollary of the above Theorem 1, which states the restriction 
on the discount factor and the corresponding achievable set for the policy based on the subset 
of MISs in Section VI. 


Corollary 1: Given the interference graph G&, if the 5 > 1 — 
vectors achieved by the policies in II approx ^ is Vapprox^ ■ 

Proof: This follows on the same lines as Theorem 1. Vamrox D 


C(G d )+1 


(Q.E.D) 
then the throughput 


= conv{7£a 


}, where 


,MIS(G d ) 


= {r(pT ),Vj G {£;(!),..., k(C(G D ) + 1)}}. It can be shown on the same lines 


that Vapprox n> is self-generating if 8 > 1 
C(G d ) + 1 extreme points. 


C(G d )+i 


. This is due to the fact that V, 


MIS(G d ) 


(Q.E.D) 


Theorem 2: For any (5 > 5 = 1 - the policy computed in Table 2 (given in Section V) 
achieves the target throughput y *(G D ). 

Proof: The policy in Table 2 (given in the paper) is based on the decomposition property 

s(Gd) Gd 

of \> mis ( g d ) (explained in the proof of Theorem 1). Define 7 (t) = a k(t)r( p 1 * ), where 

k =1 

a k (t), Wk G {1,..., s(Gd)}, V t > 0 correspond to the coefficient a. = a N+M ) at 

the beginning of time slot t in the policy in Table 2. Also, let p * correspond to the power 
vector used for transmission at time t for t > 0. First, we show that if 5 > 8 then 7 (t + 1) G 
yMis(Go), \/f > 0. Expressing 7 (t + 1 ) in terms of coefficients of a k (t), \/k G { 1 , • •■s(G'd)} as 
in Table 2. 

M* + 1) =£,*<«+1*^) = £ + — ~ s {1 ~- r(P , °°) C5) 

k =1 k=l,k^rt 

where a n (t) = nia x k a k (t). If 0 < < 1, Vk G {1, ...,s(Gd)}, k ^ r t and 0 < ° rt < 

1 then 7 (t + 1) G \? mis ( g d)' if A > max t > 0 1 — a n (t) then 7 (t + 1) G V MIS{Gd '> Vf. We know 
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that a n (t) > s(g D ) , Vt, this is true because a n (t) = max fe a k (t). The condition on S > 1 — 
implies that j(t + 1) € V MIS( ' G ' l) > 0. Therefore, we can express 

7 (f) = (1 - 5)r(p /l '‘ D ) + 8*f{t + 1 ) 

here 7 (f + 1) £ . Hence using this recursively we get, 

t 

7(0) -(!-«)£ ^r(p'° D ) = i‘ +1 7 (« + 1) 


( 6 ) 


T—0 


Here, 7 (f + 1) £ \> mis ( g d) an( j >y(o) — The above expression in (4) specifies the 

difference from the target vector and the vector achieved till time t. Next, we find the time 
T after which the norm of difference in [6] is below e. 7 (f + 1) can be bounded above (since 
yM/s(G D )j s closed and bounded), || 7 (f+l)|| < 9 bd , Vf > 0. Hence, T = - — 1 is sufficiently 

high to ensure the norm of difference in [6] is below e. Hence as e is chosen small, the policy 
would converge to the target payoff. 

Next, we give a corollary, which states the condition on the policy based on the subset of the 
MISs computed in Subsection VI-A. 

(Q.E.D) 

Corollary 2: For any 5 > 1 — c ^q d - )+1 the policy computed in Table 2 based on the subset 
of MISs computed in Subsection VI-A achieves the target throughput y* approx (G d)) 

Proof: The policy in Table 2 using the target weights cx* approx {G d) is based on the decom¬ 
position property of Vapprox ■ Here also using the decomposition property repeatedly it can be 
shown that the difference between the target and the target vector achieved till time t decreases 
exponentially which proves convergence. 

Next, we state two lemmas which will be used to prove Theorem 3 and 4. Also, note that 
any inequality between two vectors would represent a component-wise inequality. 

Lemma 1: In a bounded degree graph with N + M vertices and bound on the maximum 


degree p there exists a maximal independent set with size 


N+M 

P+1 


Proof: In a bounded degree graph the maximum chromatic number for the conventional vertex 
coloring is p + 1 ED- Let the independent sets corresponding to each color class obtained by 
using the minimum coloring with atmost p + 1 colors be given as I p+1 } and let the 

sets of the sizes corresponding to each of these independent sets be given as { 711 ,...., n p+1 }. 
Let the set of maximum size be and the corresponding size be given as rqq. We know that 
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Yli=i n i = N + M and > n t , \/i G ,p + 1}. From these conditions we can obtain 
(p + l)ri[i] > N + M. Hence, ri[ i] > 

(Q.E.D) 

Lemma 2: If the interference graph exhibits e-WNI then the expression r-fp) (as defined in 
Subsection V-C in the paper) is an e approximation to the actual rate r, (p) and this holds for 
all i G (1, N + M } and V p G V. 

Proof: To show that r -(p) is an e approximation, it is sufficient to show that max pg p |e*(p)| < 
e, where ej(p) = r-(p) — 7y(p). Solving for argmax pe p |e*(p) |, we get p t = p- nax , p 3 = 
0, Vj G A fi{Go) and p k = p™ ax , \/k G (A fi(Go) U (i}) c . This is explained as follows. 
Since the accumulative interference from non-neighbors is not accounted for in r\ (p) this means 

^_|_ djjPj _ 

r-(p) > n(p) => e*(p) > 0,Vp G V. Note that e*(p) = log 2 ( ) is 

an increasing function in p t and p k , Mh G (A fi{Go) U {i}) c and a decreasing function of 
Pj, Vj G J\fi(G D ). From this we get that e 4 (p) takes its maximum value when p, = p™ ax 
and p k = p' k ax , V/e G ( J\fi(G D ) U (i}) c and pj = 0, Vj G J\fi(G D ). The corresponding 

i+3 

= log 2 ( 


maximum value is e? 


9HV? a ^ 


1 + 


). This combined with e-WNI, i.e. 


Svfc€(Wi(G D )U{i}) c 9kiP™ aX +° 2 


Int™ ax (G D ) = Efce(M(G D )u{i}) c dkiPj 1 ™ < (2 e - l)af gives that e™ aa: < e. The same argument 
holds for any i. 

(Q.E.D) 


Theorem 3: If the constructed interference graph Go* exhibits SLI and e-WNI, then the 
proposed policy computed through Steps 1-5 of Section V-A achieves within e of the optimal 
network performance W*. 

Proof: With constraint tolerance of e, the optimization problem in Step 5 of the design 
framework is stated as follows. 
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maxW(yi(G D ), ....,y N+M (G D )) (7) 

y,a 

Vi{G D ) > R? in — e, Vi e {1 ,N + M} 

s(G d ) 

Vi(G D )= ^2 a k ri{ p 7fcD ), Vi £ {1 ,N + M} 

k= 1 

> 0, V/c € {1, ...s(Gd)} 

s(G d ) 

J2 ak = 1 

k =1 

The optimal argument of the above problem is y *{Gd)- Define y *(Gd*) = arg max D < D «s t W (y*(Gn)). 

Consider the design problem in in Subsection IV-A (given in the paper). The feasible region 
for average throughput vectors for the design problem is a subset of conv{7Z const j. This is 
explained as follows. The throughput vector achieved by an interference management policy 
can be written as: R(7r) = (1 — 5) <^ r (P*)’ h ere Rl 77 ) = (-Ri(tt), ..., Rjv+m(^))- The 

coefficients of r(p*) are positive and the sum of these coefficients is 1 and r fp') £ 1Z co " M , 
which implies that R(7t) g conv{TZ const }. Let v* £ com/{lZ const } be the optimal solution to the 
design problem in Subsection IV-A (given in the paper) with weighted sum throughput as the 
objective and the corresponding optimal value of the objective is W* = where 

Wi > 0 Vi e {1, ...,N + M} and Y^h=\ M w i = 1- Note that we have assumed that the design 
problem in Subsection IV-A is feasible, otherwise if it is not feasible then clearly the proposed 
framework will be infeasible as well. Expressing v* in terms of throughput vectors in 'JZ consl as 
follows, v* = Y2j =i where £ 'JZ const , Q- > 0, Vj £ (1,..., N + M} and 0 3 = 1. 

Here, v* > R mm where R mm = [R™ in ,R r f!i l + M \ and the inequality between the vectors is a 
component-wise inequality. Our aim is to show that there exists v £ V mis ( g d*) which is e close to 
the optimal and satisfies the minimum throughput constraint within a tolerance of e. Let r(p J ’*) = 
r 7 ’* and let r(p J ’*) £ TZ' a onst be the corresponding throughput taking only the interference from 
neighbors into account (given in Subsection V-C in the paper). Let v £ conv{1Z c ° nst } defined as 
follows v = Y2=\ @j r ' ( P' 7 ’*)- Since r is computed only from the interference contribution of the 
neighbors we have r (p J> ) > r(p J ‘*) and from SLI we know that there exists v J £ 
which satisfies 

v' J >r'(p J ’*) > r(p J ’*), Vj £ (8) 
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Let v* = J2j=i Using above ([ 8 ]) we get, 

v' > V > V* 

Expressing in terms of the throughput vectors in j ^ r7S(C '}*) we g et Y ',j — P 7fc ° 

and the corresponding actual throughput is given as \ :J = Ylk=\ } * P{ r (V lk )• From the condi¬ 
tion in the Theorem we know Gp- exhibits e-WNI which means Irit r l ,ax (Gp-) < (2 e — l)of \/i G 
{1,.., N + M}. Hence using Lemma 2, we have r-(p) is an e approximation to rj(p) and this 
holds for all i G {1, N + M}. Using r(pG ° ) > r (pG D ) — e, Vfc G {1,s(Gc*)} we get 

yJ > y'-i _ e . (10) 

Hence, using the lower bound on v J > r(p J ’*) we get v J > r(p J ’*) — e. Also, the same can be 
done in general Vj G Using this we get Yl]=\ )T^ =1 ^-/^(pG " ) > 

v* — e. Let v = Yl]=i ®j yj anc l Id v = [v i. ,.,vn+m}- We can get the following relationship. 

v > v* — e (11) 

It can be seen that we can write v = a k r (P Ik ° ) with a k = Yl]=\ Since 

C(k > 0,Vfc G {1, s(Gd*)} and Ylk=i D *^ a k = 1, which means v G \P mis ( g d*)_ This gives 
Yln=\ h ' w i v i > X/J=t A/ w i v i ~ e ant l v > R mm — e. Hence, v is a feasible throughput vector 
for 0 . Since y*(G p*) is the optimal solution to the above problem in (j7]), we can state the 
following, 

N+M N+M N+M 

Y w iVi {Gd* ) > Y WiVi - Y WiV i - e ' ( 1 2) 

i=l i= 1 z=l 

This proves the result. 

(Q.E.D) 

Theorem 4: For any homogeneous network with interference graph that has a maximum degree 
no larger than p and exhibits C-WNI with (p + 1) < min{^^^^, (1+ ^ )(c) log 2 ( 
the policy in Subsection VI-A achieves a performance W{y* approx {G D )), with a guarantee that 
W(y:„o*(G D )) > ( 1 -jWp> w^(y-(G D ), where 7 = (3(p + 1)),-=^. 

Proof: We make the following homogeneity assumption p™ ax = p rnax . R." tn = R min : ma x(D UE ) 
A. Wi = N } M - These quantities are fixed to understand the effect of scaling of network size, i.e. 
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N + M. Let r"'°' K = r, (p, = p rnax . p_ i = 0) The optimization in Step 4 in Section V-A is stated 
with weighted sum throughput as the objective. 


N+M 


max 
y ,£* 

E 

2 — 2. 

WiVi 

(13) 



s{G D ) 


subject to 

Vi = 

T+P 1 '' °), Vi G (1, ..., N + Mj 




3=1 



Vi > R min 

s(G d ) 

Y a i = 1 

3 = 1 

a j > 0, Vj G {1 ,s(Gd)} 


We formulate another optimization problem whose solution is an upper bound to the above. 

(14) 


max 
y ,<* 


subject to 


N+M 

Y w m 

i =1 

s(G D ) 

Vi= J2 

3 =1 


s(G n ) 

J2 a i = 1 

3 = 1 


a jr r%^ T ou, Vi e {!,...,iV + M} 


ie/. 


aj > 0, Vj G {1,..., s(Gd)} 


Here 1 r G„ is an indicator function which takes value 1 if i G I 


Go 
3 


and 0 otherwise. Note that 
the solving the above optimization problem in ( fT4| ) is equivalent to finding the maximum weighted 
independent set [40] with the weights assigned to each vertex i given as w r = Let 


N+M' i 
N+M 


y *(Gd) be the optimal solution to ( [13] ) and the corresponding optimal value is Wiy*(GD)- 

Let the maximum weighted independent set which is a solution to ( [T4| ) be denoted as I^ D and 
hence the optimal value of the objective in ( [T4| ) is Ejll W i r T aX1 iGl ^n- The solution to the 
problem in (jT4]) is an upper bound to the solution of (|~i~3T) , this is formally stated as 


N 

E 

i —1 


N+M 

> E 

1=1 


Wiy*{G D ) 


(15) 


This is explained as follows. Let the feasible sets of problem ( [13] ) and (JT4]) be given as T\. To 
respectively. If [cK , ,y ] G T\ then for the same a. the corresponding y" = Yl s j=l oi j r max \. 

// / D . 

satisfies y > y since r” 10,21 > ri(pL' ). 


t G D 

3 
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Next, we use the fact that if the interference graph exhibits (j-WNI, i.e. {Irrt r - iax (G [>) = 


< (2 C - l)cr 2 } then rg(p)' = log 2 (l + 


9iiV 


E T\ 

iWd# g i iP - ^ ^ J lllW11 ' ~ 1W62Vi ^ e n i{ g d) 9+pG° 2 

of rj(p). This follows from Lemma 2. Thus, we have 

rG Ds 


ien 


n (P ) < C 


is an (' approximation 


(16) 


The weight of approximate weighted maximum independent set I G ° 


E N+M m axi jG d 

i=l 1 * !)■ 1 KC[G D )+l) 

algorithm in [40], hence, we can write 

„max-t 


k(C(G D )+ 1) 18 § 1Ven aS 
is r] approximate independent set computed using the 


N+M 


i—1 


WiT, 


ieir 


> 


N+M 




WiT , 


ieC+ 


(17) 


HC(g d )+ i ) 

If «fc(C(G D )+ 1 ) = 1 and rij = 0, Vj ^ k(C(Go) + 1) then the resulting value of the objective 
function is w i r i(v k(c{GD)+1) )- From ( fT6] ) and ( [T7| ) we have, 

1 


N+M 


w i r i ('p Ik< ~ c( - G D)+ o) > 


i =1 


N+M 

E 

Z— 1 




C-l 

Ij g / G c 

* fcJ fc(C(G D ) + l) 


c> 


N+M 


l + r >' + 


(>,<-i ia «,K OS) 


The minimum value that the expression (jU- Wi+ ax 


l._ r a D ) can take is given as 




log 2 (l- 


i p" 


). To derive this, first substitute the value of Wi = 


* N+M 

N+M 


1 . From Lemma 

Also, using 


(pH-1) ( 77 -hl) AW 62V- L 1 /\jip q -2 

1 we know that there exists a maximal independent set with size no less than p+l 
the minimum value of the direct channel gain 9a > aW we get r™ ax > log 2 (l + ^^). 
Combining the fact that maximal independent set has a size no less than and r"' ax > 

. -1 „max . . . . 

log 2 (l + we get the minimum value of the expression. 

From the condition in the Theorem we have ( < ( - p+1 ^ T?+1) log 2 (l + where 0 < k < 1 

determines the distance from the optimal solution. If ( is selected based on this threshold then, 

Yld=i M w i r i{ p y '= (c(G D )+1) > Y+iN^=\ M w i r T ax ^ . gI g d)- Using the fact that solution to ( [T4l ) is 


an upper bound to ( fT3] ) as stated in ( fi3| ), we get the following. 


N+M 

E 

i =1 


WjA(p /fe(c(Gl3)+1) ) > 


N+M 

K . \ - 


1 + r i ti 


uiji > 


N+M 

K \ -« 


1 + 7) ^ 


WiV* (G 


D 


(19) 


1=1 


Next, we state the optimization problem to compute W (y* approx (G d)) which is similar to the 
problem in Step-4 in the Section V-A but uses the subset of the MISs computed in Subsection 
VI-A, here y * approx {GD) is the corresponding optimal argument. Note that W (y* approx (G d)) is 
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the optimal value that can be attained by the policy based on the subset of MISs computed in 
Subsection VI-A in the paper. 

N+M 

max V Wiyi (20) 

y ,ol zJ 

i= 1 

c(g d )+i Gd 

subject to i/i = E ajTii pV)), Mi G {1,...,N + M} 

3 = 1 

Vi > R min , Mi G {1, N + M} 

C(G D )+1 

I] «i = 1 

J =1 

a j > 0, Vj G {1, C(Gd) + 1} 

Using r/ M > and the fact that Int™ ax (G d) < (2^ — l)cr 2 (£-WNI) ensures that following 


is true 


r ( (p r ‘M) > log 2 (l + P 


A „p 2 c ( 7 -j) 1 i 6 /g;- v J e < 1 ....,C(G D )} 

We now develop a feasible solution to the above problem in ( |20| ). Assign 3, 


( 21 ) 


approx,k(j) 


R n 


Vj G {1, C(Gd)}- Hence, we can write 


( 22 ) 


l0g 2( 1+ A"P2C CT 2) ' 

C(G d ) g C(G d ) 

Y Approx, kuMV Ik &)>R min Y l ieI G D ,Vje{l,...,C(G D )} 

3=1 ' 3=1 kU) 

The maximal independent sets {/^,4 ^c(g d ))} are obtained by adding vertices in a greedy 
manner in ii) to the independent sets constituting color classes, as in Subsection VI-A. Hence, 


all these MISs together contain all the vertices in the graph, this ensures J +7 


c(g d ) 


3=1 


~ L 


Therefore, this assignment of f3 apP rox,k(j) ensures the minimum throughput constraint in ([20]) is 


satisfied. However, the it still needs to be checked if /3 apP rox,k(i) < 1. Since C(Go) < 


fha 


3.(p +1) from [38]. The condition on p, i.e. p < — 1 ensures Yl < ji < i D ' > A 


K" 


j=i r' approx ,k(i) 

< 1. Hence, the solution obtained is feasible. The remaining fraction 


V (3(p+l)); . \ 

10g2(.l+ A np 2eCT 2) 

1 — 7 is assigned in such a way to ensure a constant competitive ratio. 


Assign 3a Ppr0 x,k(C(G D )+ 1 ) = ( 1 - 7 ) then the lower bound is given as W(y* pprox (G D )) > 
(1 — 7 ) Yl!i=i M 'iCiU(p 7fc(C(G£>:)+1) )- This combined with the lower bound in ( fT9[ ) derived above 
gives: 


(1 - «)(1 - 7 ) 


N+M 


W(y‘ appr „(Go)) > " V WivKGo) (23) 


1 + 7 


i =1 
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(Q.E.D) 


Discussion on Coloring based policies vs MIS based policies 

For this discussion we assume non-neighboring UEs in the interference graph do not interfere 
at all. In a coloring based policy, the set of UEs that are scheduled form an IS of the interference 
graph, which need not be maximal. Each such IS which is not maximal can be extended to 
MIS. The throughput vector obtained by scheduling UEs in a particular IS is dominated by the 
throughput vector corresponding to the MIS obtained by extending the same IS. This is because 
the non-neighboring UEs which are added in order to extend IS to MIS will now have a positive 
throughput and the UEs already in the IS will not be affected because the non-neighboring UEs 
do not interfere. 
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